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This  paper  deals  with  a  comparative  study  of  circle  criterion  based  nonlinear 
observer  and  Hrj  observer  for  induction  motor  (1M)  drive.  The  advantage  of 
the  circle  criterion  approach  for  nonlinear  observer  design  is  that  it  directly 
handles  the  nonlinearities  of  the  system  with  less  restriction  conditions  in 
contrast  of  the  other  methods  which  attempt  to  eliminate  them.  However  the 
//*>  observer  guaranteed  the  stability  taking  into  account  disturbance  and 
noise  attenuation.  Linear  matrix  inequality  (LMI)  optimization  approach  is 
used  to  compute  the  gains  matrices  for  the  two  observers.  The  simulation 
results  show  the  superiority  of  Ha  observer  in  the  sense  that  it  can  achieve 
convergence  to  the  true  state,  despite  the  nonlinearity  of  model  and  the 
presence  of  disturbance. 
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1.  INTRODUCTION 

For  a  long  time,  the  induction  motor  is  considered  as  a  principal  workhorse  in  the  industry  due  to  its 
robustness,  high  reliability,  relatively  low  cost,  modest  maintenance  requirements  and  efficiency  [1,  2]. 
However,  induction  motor  is  also  known  as  a  complex  nonlinear  system,  in  which  time-varying  parameters 
entail  additional  difficulties  for  machine  control,  conditions  monitoring  and  fault  diagnostic  purposes  [1], 
The  main  problem  with  induction  motor  industrial  applications  is  that  only  a  few  state  variables  of  the 
machine  are  available  for  on-line  measurements.  This  is  due  to  technical  and/or  economical  constraints. 

In  order  to  perform  advanced  control  techniques,  conditions  monitoring  and  faults  diagnosis  there  is 
a  great  need  of  a  reliable  and  accurate  estimation  of  the  key  unmeasurable  state  variables  of  the  machine.  In 
this  context,  the  observer  design  theory  seems  to  be  an  ideal  solution.  Over  the  last  two  decades,  nonlinear 
observer  design  problem  has  received  much  attention  in  the  literature.  Several  attempts  have  been  made  for 
particular  classes  of  nonlinear  systems.  Existing  approaches  can  be  roughly  classified  as  follows:  Nonlinear 
state  linearization  approaches,  high-gain  observers,  geometric  algorithms,  variable  structure  design 
procedures,  and  algebraic  techniques  [3-6]. 

The  problem  with  sensorless  induction  motor  industrial  applications  is  which  approach  that  provide 
the  most  accurate  and  reliable  estimation  of  unmeasurable  state  variables  of  the  machine  system. 

In  this  paper  we  focus  our  attention  on  the  performance  comparison  of  the  so  called  circle  criterion 
approach  and  H,rj  observer  design  for  induction  motor  (IM)  system.  The  advantage  of  the  circle  criterion 
approach  is  the  direct  handling  the  system  nonlinearities  by  exploiting  their  properties.  This  approach  is  less 
restrictive  compared  to  the  other  methods  as  linearization  approach,  high  gain  observer  and  variable  structure 
design  procedures  witch  attempt  to  eliminate  them  [5,  6].  The  H«,  observer  is  used  to  guaranties  the 
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robustness  against  nonlinearities  assumed  to  be  equivalent  to  some  uncertainties  due  to  internal  and  extremal 
disturbances  and  measurement  noises  [7-12]. 

To  make  a  comparative  study  between  the  circle  criterion  based  observer  and  observer  design  for 
induction  motor  derive,  we  use  one  of  the  standard  model  of  an  induction  motor.  This  type  of  nonlinear 
model  is  generally  used  for  perfonning  nonlinear  control  strategies,  conditions  monitoring  and  faults 
diagnosis  of  electric  induction  machine  systems. 

The  paper  is  organized  as  follows:  In  the  second  and  third  section  we  present  the  theory  of  nonlinear 
observer  and  the  nonlinear  observer  based  circle  criterion  respectively.  The  robust  H „  observer  is  presented 
in  the  fourth  section.  In  the  fifth  section  we  present  the  considered  nonlinear  induction  motor  model.  Finlay, 
we  present  simulation  results  and  comments.  A  conclusion  ends  the  paper. 


2.  THE  NONLINEAR  LUENBERGER  OBSERVER  DESIGN 

We  recall  that  an  observer  is  a  dynamical  system  which  uses  the  available  input-output  data  to 
reconstruct  the  unmeasurable  system  state  variables.  It  is  a  “soft  sensor”  that  plays  an  important  role  not  only 
in  sensorless  control  techniques  but  also  in  conditions  monitoring,  fault  diagnosis,  predictive  maintenance 
and  fault  tolerant  control  techniques  [1,2]. 

In  this  paper,  we  consider  a  class  of  nonlinear  systems  that  can  be  decomposed  in  linear  and 
nonlinear  parts  as: 

x(t )  =  Ax(t )  +  F  (x(t),u(t))  ^ 

y(t)  =  Cx{t)  m 


Where  stand  for  the  system  state  variable,  U ^  for  system  input  and  ^ ^  for  the  system  output. 

The  function  represents  the  system  nonlinearities.  For  this  type  of  nonlinear  system,  a  general 

nonlinear  observer  (Luenberger)  expression  is  as  follows: 


x(t)  =  Ax(t)  +  F(x,  u )  +  L(y(t)  -  y(t)) 
y(t )  =  Cx(t) 


Where  x^  stands  for  the  estimated  state.  We  assume  that  the  nonlinear  function  F(-X(t),u(t))  js 
locally  Lipschitz  with  respect  to  the  state  variable  x^  . 

II  With  ^  >  ^  as  the  Lipschitz  constant; 


3y/||A(x(0,t<t))-77(x(040|^y|K0-^(0||Wl-th  /  >  0, 


Theorem  1  [13]:  Consider  Lipschitz  nonlinear  system  (l)-(2)  along  with  the  observer  (3)-(4).  The  observer 
error  dynamics  are  asymptotically  stable  with  maximum  admissible  Lipschitz  constant  if  there  exist  scalers 

s  >  0  and  b ’  >  and  matrices  ^  >  ®  and  ^  such  that  the  following  LMI  optimization  problem  has  a  solution 
min(  £) 


(A-LC)t  P  +  P(A-LC) 
P 


P 

P 


<  0 


2  y 


(5) 


1 

2 


4 


p 


With  r*=max(r)  =  ^ 


(6) 
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Proof:  Defining  the  observer  error  as  e ^  x^ 

e(t )  =  (A-  LC)e  +  F(x(t),  u(t))  -  F(x(t),  u(t)) 

In  the  following  we  use  notations: 

F  =  F(x(t),u(t ))  and  F  =  F(x(t),u(0) 

The  time  derivative  of  the  Lyapunov  function  1  = 
V(t)  =  eT  [(A  -  LC)T  P  +  P(A  -  LQ]e  +  2  eT  P(F  -  F) 
By  using  the  Lipschitz  property  we  have: 

=>  2||e7>(F  -  F)||  <  2y\ef  P||e|| 

Then  equation  (8)  becomes: 

V(t)  =  er[(A  -LC)tP  +  P(A-LQ  +  2  y;rP\e  <  0 
Those  leads  to: 


then  the  observer  error  dynamics  are  given  by: 

(7) 

Pe  can  be  expressed  as: 

(8) 

(9) 

(10) 


(. A-LQTP+P(A-LQ  +  2yPp-lP<0 


(ID 


By  using  the  Schur  complement  theorem  we  therefore  obtain  the  following  linear  matrix  inequality  (LMI): 


(A-LC)T  P  +  P(A-LC)  P 

P  __  J— 

2r  _ 

The  Lipschitz  constant  can  be  computed  as  follows: 
For  any  Lyapunov  function  F  =  e  Pe  we  dave: 

W^IHI2  ^eTPe<Xm^{P)\ef 


(12) 


(13) 


2reTPe<2yAmin(P)\\e\\- 

Let:  (A  ~LC)Tp  +  p(A  ~  LC)  =  ~Q 

Then  the  time  derivative  can  be  written  as: 

V(t)  =  -eT  Qe  +  2yP  ||e||2  <  0 

F(0  =  -^dn(G)|Hf  +  2yXmax(P)\\ef  <  0 

So:  2r^maJp)  <  Amm  (Q) 

2Tmax(P) 

In  addition,  since  P  is  positive  definite,  suppose 


(14) 

(15) 

(16) 

(17) 

(18) 

I  .From  the  equation  (18)  we  can  write: 
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\-2yXmm{P)>Q 

And  a{P)  =  {P) 

So  from  (19): 

2/ 

Which  is  equivalent  to: 

(—  )2I-PTP>0 
2  r 


Using  Schur  complement  lemma 


2y 


I 


P 


Defining  Y  ,  (6)  is  achieved. 


(19) 


(20) 


(21) 


(22) 


3.  CIRCLE  CRITERION  BASED  NONLINEAR  OBSERVER  DESIGN 

In  contrast  of  the  linearization-based  and  high-gain  approaches  which  attempt  to  eliminate  the 
system  nonlinearities  using  a  nonlinear  state  transformation  or  to  dominate  them  by  a  high  gain  term  of 
correction,  circle-criterion  exploits  the  properties  of  the  system  nonlinearities.  In  its  basic  form,  introduced 
by  Arcak  and  Kokotovic  [14],  the  approach  is  applicable  to  a  class  of  systems  that  can  be  decomposed  in 
linear  and  nonlinear  parts  with  a  condition  that  the  nonlinearities  satisfy  the  sector  property  [14-16]. 


3.1.  Basic  sector  properties 

A  memoryless  nonlinear  function  F^z’  ^ ^  +  ^  [x^  R  is  said  to  belong  to  the  sector  ^  +00[ 
if  0  Let  Vl  and  Vl  two  real  positive  numbers,  by  setting  z~vi  v2  ancj 

F(~P)  -  [F(vut)  F(y2,t)\ ,  ^  a^ove  sector  property  is  equivalent  to: 


(V[  -  V2)[F(v1;0  -F(v2,0]  -  0  VV|,V2  e  R+ 


(23) 


Relation  (23)  states  that  the  nonlinear  function  ^z’  ^  is  a  nondecreasing  function.  On  the  other 
hand  if  a  continuously  differentiable  function  the  above  relation  is  equivalent  to  [7],  [8]: 

—F(z,t)>  0  Vzeli 

dz  (24) 


F(z,t) 


If  the  nonlinear  function 
such  that: 


F(z,t) 


does  not  satisfy  the  positivity  condition  (24)  we  introduce  a  function 


F(z,  t )  =  F(z,  t)  +  pz. 


P >■ 


4~F(z,t) 


,VzeU 


(25) 


And: 
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—F(z,t)  =  —F(z,t)  +  p>  0  Vz  < 
dz  dz 


(26) 


In  the  multivariable  case  the  sector  property  can  be  written  as:  z  F(zf)  >  ^  .  Where  z  and^z’^ 
are  respectively  vectors  of  an  appropriate  dimension. 

3.2.  Nonlinear  Observer  Design 

The  circle  criterion  based  nonlinear  observer  design  can  be  performed  for  a  class  of  nonlinear 
system  that  the  model  can  be  decomposed  into  linear  part  and  nonlinear  part  as  the  following  [14-16]: 


x(t )  =  Ax(t)  +  </[u(t),  y(t)\  +  GF[Fl.x(t)\ 
y(t)  =  Cx(t) 


(27) 

(28) 


Where  ^  ^  and  ^  are  known  constant  matrices  with  appropriate  dimensions.  The  pair  is 

assumed  to  be  observable.  The  term  74>'!f),w(t)]  js  an  ari-,itrary  real-valued  vector  that  depends  only  on  the 
system  measured  control  inputs  and  outputs  .  The  nonlinear  part  of  the  system  is  included  in  the 

term  which  is  a  time -varying  vector  function  verifying  the  sector  property.  In  the  following  we 

recall  the  main  theorem  and  conditions  that  are  used  in  this  work  to  study  the  feasibility  of  nonlinear  observer 
design  for  induction  motor  sensorless  control  with  respect  of  circle  criterion  or  sector  property.  A  detailed 
proof  of  the  theorem  is  presented  in  reference  [17]. 

Theorem  2  [14,  15]:  Consider  a  nonlinear  system  of  the  form  (27)-(28)  with  the  nonlinear  part  satisfying  the 


circle  criterion  relations  (23)-(26).  If  there  exist  a  symmetric  and  positive  definite  matrix 
Kg  R.p 

of  row  vectors  such  that  the  following  linear  matrix  inequalities  (LMI)  hold: 

(A-  LC)T  P  +  P(A-  LC)  +  Q  <  0 


P  eRn 


and  a  set 


(29) 

(30) 


PG  +  (H-KCy  =0 

The  nonlinear  observer  design  refers  to  the  selection  of  the  gain  matrices  L  and  K  satisfying  the  LMI 

conditions  (29)-(30).  With  ^  ~  1 

£  is  a  small  positive  real  number. 

Then  a  nonlinear  observer  can  be  designed  as: 


conditions  (29)-(30).  With  @  £^n  as  a  defined  positive  known  matrix,  is  an  n-th  order  unity  matrix  and 


x(t)  =  Ax(t)  +  <f{ii(t),y(t)] + Hy(t)  -  y(t)\ + GF[Hx(t) + K(y(i)  -  y(t))] 


y(t)  =  Cx(t) 


(31) 


(32) 


And  e(0  -  *(0  x(0  -  ,  where  x^is  the  estimate  of  the  state  vector  of  the 

nonlinear  system.  One  can  see  that  the  structure  of  the  nonlinear  observer  is  composed  of  a  linear  part  that  is 
similar  to  linear  Luenberger  observer  and  a  nonlinear  part  that  is  an  additional  term  that  represents  the  time- 
varying  nonlinearities  satisfying  the  sector  property. 


Proof:  The  dynamics  of  the  state  estimation  error 


e(t )  =x(t )  -  x(t) 


are  given  by: 


e(t)  =  (A-  LC)e(t )  +  G[  f(H.x(t))  -  F{H.i(t )  +  K(y(t)  -  y(t))} 


(33) 
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Let  v,  =  H.x(t)  ,mdV2  =  H:X(t)  +  K(y(t)-y(t))  ^  by  scttingz  = vi  “v2  =  ■ ~KC)e(t)  >  the  tenn  between 

brackets  in  (33)  can  be  seen  as  a  function  of  the  variable  z  then:  ^Vl  ^  2  ^  -  7  ( z ’  ^  . 

The  expression  ^Vl  -  )  -  F(v2)\  -  -F(z,t)  satjsges  tbe  pr0perty  0f  the  sector. 

Taking  into  account  the  above  result,  the  error  dynamics  in  (33)  can  be  rewritten  as: 


e(t)  =  (A-  LC)e(t)  +  G.F(z,  t ) 


(34) 


z  =  (H-KC)e(t ) 


(35) 


Relations  (34)-(35)  show,  once  again,  that  the  error  dynamics  can  then  be  considered  as  a  linear 

system  controlled  by  a  time-varying  nonlinearity  function  satisfying  the  sector  property.  Circle 

criterion  establishes  that  the  feedback  interconnection  of  a  linear  system  and  a  time -varying  nonlinearity 
satisfying  the  sector  property  is  globally  uniformly  asymptotically  stable  [14,  15], 

Based  upon  the  error  dynamics,  relation  (34)-(35),  the  nonlinear  observer  design  problem  is  then 
equivalent  to  stabilization  of  the  error  dynamics  problem.  To  this  end  a  candidate  Lyapunov  function 

V  =  e  Pe  js  considered.  In  order  to  ensure  asymptotic  stability  of  the  observer,  the  derivative  of  the 
candidate  Lyapunov  function  must  be  negative.  With  the  help  of  relation  (34)  and  (35)  the  derivative  of  the 
Lyapunov  function  becomes: 


V  =  er[(A-LC)rP  +  P(A-  LC)\e  +  FT(z,  t)Gr  Pe  +  erPGF(z,t) 


By  setting: 


(A-LCY  P  +  P(A  -  LC)  <  -Q 


(36) 


(37) 


And 


PG  =  -(H-KCy 


(38) 


With  @  ^ n  and  £  >  ®  ,  the  derivative  of  the  Lyapunov  function  can  be  rewritten  as: 

V<-eTQe-2.zT  F(z,t)  (39) 

Thus  ends  the  proof. 

Note  that  the  existence  of  observer  (31)-(32)  is  conditioned  by  the  solution  of  LMI  conditions  (29)-(30).  By 
solving  LMI  constraints,  observer  gain  matrices  L  and  K  that  guarantee  observer  convergence  are  then 
computed.  In  Ibrir  [16],  the  author  has  investigated  the  study  of  globally  Lipschitz  systems  and  bounded-state 
nonlinear  systems.  Bounded-state  nonlinear  systems  constitute  a  large  class  of  system  that  includes  electric 
machine  systems.  Electric  machine  models  involve  the  magnetic  flux  as  a  key  and  bounded  state  variable  that 
combined  with  other  state  variable  of  the  machine,  such  as  rotor  angular  velocity,  leads  to  the  nonlinear  part 
of  the  machine  model.  This  is  due  to  the  effect  of  the  magnetic  material  saturation  property  that  is  similar  to 
the  sector  nonlinearity. 


4.  NONLINEAR  H,  OBSERVER  SYNTHESIS 

In  this  section,  we  propose  to  study  the  robust  Hx  observer.  This  observer  is  used  in  the  case  of 
systems  involving  uncertainties  in  the  model  measurements  [7,  9,  10,  1 1,  12,  18]. 

A  literature  review  shows  that  many  useful  H ^  filtering  approaches  have  been  developed  for  several  kinds  of 
systems.  One  can  cite  the  Hro  observers  for  one-sided  Lipschitz  nonlinear  systems,  the  Ha 0  observer  for 
singular  Lipschitz  nonlinear  systems  and  the  Hro  control  and  filtering  for  uncertain  Markovian  jump  systems 
with  time -varying  delay  [9,  10,  11,  19]. 

The  aim  of  the  Hx,  approach  is  to  design  a  full-order  filter  such  that  the  corresponding  filtering  error 
system  is  asymptotically  stable  and  satisfies  a  prescribed  Hx  level  disturbance  attenuation.  By  using  a 
Lyapunov  function,  sufficient  conditions  are  formulated  in  terms  of  linear  matrix  inequalities  (LMIs). 
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To  do  this,  we  extend  the  result  of  the  previous  section  to  nonlinear  robust  Ho o  observer  design  method.  We 
consider  the  system  (l)-(2)  with  an  additional  disturbance  term  as  follows: 

x(t)  =  Ax(t )  +  F(x{t),  u  (t))  +  Bd  w(t )  (40) 

y(t)  =  Cx(t)  . 


Where  vi  e^2^’00)  is  an  unknown  exogenous  disturbance. 

x(t') 

The  goal  is  to  rebuild  the  state  v  ’  of  the  system  (40)-(41)  with  some  accuracy  despite  the  presence  of 
disturbance  term  "  ^  . 

We  suppose  that:  z  _  ^eit) ,  wjtjj  H  a  known  constant  matrix  such  that  ^  ;  with 

The  corresponding  observer  for  the  system  (40)-(4 1 )  is  given  by: 


x(t)  =  Ax(t)  +  F(x(t),u(t ))  +  L[y{t)  -  J>(f)] 


(42) 


y{t)  =  Cx(t ) 


(43) 


Our  purpose  is  to  design  the  observer  parameter  L  such  that  the  observer  error  dynamics  are  asymptotically 

JJ 

stable  and  the  following  specified  norm  00  upper  bound  is  simultaneously  guaranteed. 


Lemma  1:  For  any S  ^  and  any 
2  xTy<xTPx  +  yTP~ly 


positive  definite 


matrix 


P  e  R,,xn 

9 


we  have: 


Theorem  3  [13]:  Consider  stochastic  Lipschitz  nonlinear  system  (40)-(41),  and  the  corresponding  observer 
(42)-(43).  The  observer  error  dynamics  are  asymptotically  stable  with  the  minimum  of  the  norm  and  ^  ;  if 

there  exist  scalers  01  s>  1,  £  >  ®  and  h  >  ®  and  matrices  ^>>®and  ^  such  that  the  following  LMI 
optimization  problem  has  a  solution. 
min(  £) 


HTH  +  Uy  +  --as)I  PBd 

2  y 

btap  -z 

'1-^7  7> 

2r  <0 

p  \-F2(H) 

2r 

With  L  =  P-'G  and  =  ^ 

Proof:  The  observer  error  dynamics  are: 

ft)  =  (A-  LQe(t)  +  {F  —  F)  +  Bdw(t) 


The  time  derivative  of  the  Lyapunov  function  e  Be(t)  gjven  by: 

V(t)  =  eT[(A  -  LC)tP  +  P{A  -  LC)]e 
+  2 erP(F  -F)  +  erPBdw(t)  +  w(t)T BT,Pe 


(44) 


(45) 


(46) 


(47) 
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Using  Lemma  1  and  the  Rayleigh  leads  to  the  following: 

2 erP(F  -F)<  ePe +{F-  F)T  PP~lP(F  -  F) 

(48) 

2eT  P(F  -  F)<  eTPe  +  (F- F)T P(F  -  F) 

(49) 

The  right  term  of  (49)  can  be  written  as: 

eTPe<Xm,x(P)\\ef  =XxmJP)eTe 

(50) 

\\(F-F)P(F-F)f  <r2Xmax(P)\\4 

(51) 

The  left  tenn  of  (49)  can  be  written  as: 

l\eT  P(F  -  F)\<2yXmax(P)\\ef 

(52) 

Replacing  (5 1  )-(52)  in  equation  (49),  leads  to: 

2r4,axCP)H2  -  4iax(-P)  H  2  +r2^max(^)  H 

Assume  ^  ^  ,  in  the  deterministic  case  we  have  w 

function  Lyapunov  becomes: 

V(t)  <  -eTQe  +  2 eT  (F  -  F)  <  0 

2 

=  0  and  the  corresponding  time  derivative  of  the 

(53) 

V(t)<-eTIe  +  2rXmm(P)\H2  <0 

(54) 

Then: 

l-2M,nax(n>0 

(55) 

a{P)  <  1- 

2  r 

(56) 

2eTP(F-F)<(l  +  r2)Xmax(P)eTe 

(57) 

2erP(F  -  F)  <  (1  +  y2)  —  (P)eTe 

2  r 

1  1  T 

^  -(/  +  ~)eTe 

2  y 

(58) 

In  the  stochastic  and  general  case  of  w  ^  ®  and  ^ 
becomes: 

V(t)  =  eT[(A  -  LC)T P  +  P{A  -  LC)]e 
+  2 eTP{F  -  F)  +  eT PBdw(t)  +  w(t)TBdPe 

the  above  equation  of  Lyapunov  derivative 

(59) 

V (0  ^  ~Ana ykQV e  +  ^  0  +  ~)eT e 

2  Y 

+  wr(t)B'dPe+erPBdMt) 

(60) 
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V(t)=  -  a seT e  +  —  (y  +  —  )eT e 

2  y 

+  wT  (t)BdPe  +  eT PBdw(t) 

(61) 

V(t)  =[—  (\+r)-a)Y e+vf  ( t)BrdPe+erPBdw(t ) 

2  y 

(62) 

z  Heft )  h  [sa  known  constant  with  HI  H  M  1 

Now,  we  define  the  following  criterion  to  minimize  the  difference  between  the  energy  of  the  estimation  error 

and  the  energy  inducted  by  the  disturbances: 

00 

J  =  J (zTz  -  fyvT w)dt 

0 

(63) 

Thus: 

00 

J  <  j"  (zT  z  -  <fwT  w  +  V ( t))dt 
0 


(64) 


So  a  sufficient  condition  for  J  <0  js  tha.t: 

\/t  e  [0,oo),  zT z- C>vT w +  V(t)  0 

This  means  that:  “  fw  w<0  ^ccause  V{t)  js  a  negative  function)  hence: 


||z||2  <  ^||wf  =>  ||z||  <  V^IHI- 

But  we  have: 

zTz  -  <ffw+  V(t)  =J H1 He-  & w+  V(t ) 

1 


<<?HTHe+[—(\+f)-a,)Ye 
2 y 

+e PIfaw+w B^Pe—fyf' w 


Then,  we  can  write: 


\  T  T 1 

HTH  +  d-(l  +  yl)-a,y  PBC, 

e 

je  w  J 

2  Y 

w 

BTdp  -a_ 

So  a  sufficient  condition  for  J  <  ®  is  that  the  following  be  negative  definite: 


HT  H  +  (1  +  y2 )  —  a  s  ]/  PBd 


2  y 


BTdP 


<0 


According  to  the  Schur’s  complement  lemma,  (67)  is  equivalent  to: 
HtH  +  [2 Mmax (P)  -as\  +  -PBdBT,P<  0 

I 


.-&<  0 


(65) 


(66) 


(67) 


(68) 


For  equation  (68)  to  be  negative,  the  sum  of  the  first  and  second  tenns  must  be  negative  because  the  third 
term  is  always  positive: 


(69) 
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But  as  for  any  other  symmetric  matrix,  for  H  H  ,  we  have: 

An  m(.HTH)<HTH<Amax{HTH) 


(70) 


*iax(H)I>HTH>cr2mm(H)I 

Or  according  to  the  definition  of  singular  values 

a2  (H)  +  2yAmax  (P)  -as<  0 


Or: 


AnaxO P)< 


as  -  a  2  ( H ) 

2  r 


This  is  equivalent  to  (45). 


(71) 

(72) 


(73) 


5.  SIMULATION  RESULTS  AND  COMMENTS 

To  perform  the  comparative  study  between  the  two  nonlinear  observers  previously  view,  we  use  the 
model  of  an  induction  motor.  Described  by  the  following  nonlinear  differential  equations  with,  the  stator 
current,  rotor  flux  and  rotor  angular  velocity  as  selected  state  variables  of  the  machine  [  1 ,  20,  21]. 

d  .  p  .  1 

~  ha  -  ~7  ha  +  —  <Pra  +  P  +  ~ 11  sa 

clt  1  y  ats 

d  .  P  1 

—  hP  =  ~Yhp  -P°}r(Pra  +  ~  <PrP  +  ~ 
at  1  r  ats 

d  m  .  1 

~7~(Pra  ~  ~ZTlsa  —  -prVra  ~  ^rVi-p 
at  ir  ir 

dm.  1 

~J^‘sd+  ~~<PrP 

——cor  —  a(.<praisp  —  <prpisa ) —  kjtOy  —  kpTj 


Where: 


Q„  =  ®L  a  =  ^  =  l[lz£:]  =  I^_  <T  =  1 

«„  Jl,.  m  l  cr  J  a  Llr 


—  k  =—  k  T  =  — 

Ur  f  J  ^  J  '  K 


The  indexes  5  and  r  refer  to  the  stator  and  the  rotor  components  respectively.  ,  ^  and  u  respectively 


A  , 

denote  the  stator  current,  the  rotor  fluxes,  the  supplied  stator  voltage,  is  the  resistance,  '  is  the 

T  T 

inductance,  m  is  the  mutual  inductance.  *  and  r  are  the  stator  and  the  rotor  time  constant  respectively. 

O')  .  f  .  j  /X  p 

r  is  the  rotor  angular  velocity,  J  r  is  the  friction  coefficient,  J  is  the  moment  of  inertia  coefficient,  1  is 

Q  T 

the  number  of  pair  poles,  r  is  the  mechanical  speed  of  the  rotor  and  finally  1  is  the  mechanical  load 

torque. 

This  type  of  nonlinear  model  is  generally  used  for  performing  nonlinear  control,  conditions 
monitoring  and  faults  diagnosis  of  electric  induction  machine  systems.  Performing  these  techniques  requires 
estimating  unmeasured  rotor  flux  linkage  and  rotor  angular  velocity  state  variables  based  on  the  stator  current 
and  voltage  measurements  [21]. 

To  carry  out  the  comparative  study  between  the  two  observers  previously  view,  we  use  the  induction 
motor  with  the  following  characteristics  details  on  Table  1. 
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Note  that  the  starting  conditions  X0=  [4  0  00  0]  for  the  initial  state  vector  of  the  machine  are  taken 
into  account  in  simulations  of  the  two  observers. 

In  the  application  of  the  nonlinear  observer  based  on  the  criterion  of  circle  approach,  we  must  write 
the  model  of  induction  motor  (74)  in  the  fonn  (27)-(28),  taking  into  account  the  properties  (29)  and  (30). 


Table  1  Characteristics  of  the  induction  motor 


Symbol 

Quantity 

Numerical  value 

p 

Power 

1.5  KW 

/ 

Supply  frequency 

50  Hz 

u 

Supply  voltage 

220  V 

np 

Number  of  pair  poles 

2 

K 

Stator  resistance 

4.850  £2 

K 

Rotor  resistance 

3.805  £2 

K 

Stator  inductance 

0.274  H 

K 

Rotor  inductance 

0.274  H 

m 

Mutual  inductance 

0.258  H 

('f 

Rotor  angular  speed 

297.25  rd/s 

J 

Inertia  coefficient 

0.03 1  ^2*  ^ s 

fr 

Fiction  coefficient 

0.001 14  N.s/rd 

T, 

Load  torque 

5  N.m 

The  nonlinearities  of  the  machine  system  are  function  of  the  flux  state  variable  that  is  a  bounded  state 
variable.  The  nonlinearities  of  the  model  are  of  the  form  G},  (^rd  that  can  be  expressed  as: 


«>/PrJ  =  ( G>r<Prd  +  P<»r)  ~  P°>r 


(75) 


One  can  verify  that: 


<P,-d  +  P^r)  =  <Prd  +  P  ^  0 

dcor 


With 


Ml-2 


,  then  one  can  choose 


P  =  2 


(76) 


The  first  step  of  the  simulation  consists  of  resolving  the  LMI  conditions,  relation  (29)-(30),  using  an 
adequate  LMI  tools  such  as  the  LMI  tool-box  of  the  Matlab  software.  The  obtained  nonlinear  observer  gain 
matrices  L  and  K  are  the  following: 


- 1.6749 

0.1188  " 

0.1188 

-1.6749 

L  = 

-0.7172 

-0.1075 

-0.1075 

-0.7172 

1.6201 

-1.6201 

6037  ] 

^3  =  [0.3948  -03 

=  [-0.9193  0.3948] 


The  corresponding  Lyapunov  matrix  for  this  LMI  feasibility  test  is: 


A  comparative  study  of  nonlinear  circle  criterion  based  observer  and  Ho o  observer (Farid  Berrezzek) 


1240  □ 


ISSN:  2088-8694 


P= 


0.1550  -0.07 1 0  00514  0.1486  0.0274 

-0.0710  0.1550  0.1486  0.0514  -0.027^ 

0.0514  0.1486  5.6010  0.4659  -0.0505 

0.1486  00514  0.4659  5.6010  0.0505 

00274  -00274-00505  0.0505  0.0173 


With:  e  =  °-04  . 

The  second  step  of  simulation  consists  of  injecting  the  obtained  numerical  values  of  the  gain  matrix 
and  the  vectors  in  an  S-function-based  Matlab  program  that  interacts  with  the  Matlab  Simulink  software  to 
simulate  the  nonlinear  system  and  the  nonlinear  observer. 

The  simulation  results  of  the  designed  nonlinear  observer  based  circle  criterion  are  presented  in  the 
following.  Figure  1  and  Figure  2  show  the  measured  and  estimated  stator  current  and  rotor  flux  components 
respectively.  The  Figure  3  and  Figure  4  show  the  measured  and  estimated  rotor  angular  velocity  and  the 
corresponding  load  torque  respectively  with  corresponding  estimation  error.  One  can  see  that  the  estimated 
state  variables  of  the  machine  follow  the  desired  trajectories. 

To  highlight  these  results  a  load  torque  is  introduced  in  the  simulation  at  time  of®-^sec,  the 
simulation  results  show  that  all  the  state  variables  of  the  machine  are  modified  accordingly.  Thus 
demonstrate  the  effectiveness  of  the  circle  criterion  based  nonlinear  observer  design  for  the  induction 
machine  system  state  estimation. 


. 

- 1  -  -  r  -‘-71-  -  t  -  -  i - 1  -  -  i - 

V  /ffi  A  ll  .  |  |  | 

jjstrxEES 

i  i  i  i  i  i  i 

- 4 - 

Figure  1.  Measured  (blue  line)  and  observed(green 
line)  d-stator  current  components  and  the 
corresponding  estimation  error 


Figure  2.  Measured  (blue  line)  and  observed  (green 
line)  r/-rotor  flux  components  and  the  corresponding 
estimation  error 


* 
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Figure  3.  Measured  (blue  line)  and  observed  (green 
line  )  rotor  angular  velocity  of  the  (IM)  and  the 
corresponding  estimation  error 


Figure  4.  Measured  (blue  line)  and  estimated  (green 
line)  electromechanical  torque  of  the 
(IM)  and  the  corresponding  estimation  error 
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In  simulation  experiments  of  nonlinear  “  observer  we  must  write  the  model  of  induction  motor 
(74)  in  the  form  (40)  -  (41). 

As  in  the  first  observer  approach,  the  simulation  of  the  nonlinear  “  approach  is  performed  in  two 
steps,  the  first  simulation  step  consists  of  resolving  the  LMI  conditions,  relation  (44)-(45),  using  an  adequate 
LMI  tools  such  as  the  LMI  tool-box  of  the  Matlab  software. 

To  ensure  robustness  against  nonlinear  uncertainty,  we  use  the  theorem  1  to  maximize  the  admissible 
Lipschitz  constant  Y  and  then  the  theorem  3,  to  minimize  ^  for  the  maximized  ^  . 
y*  =  0,7  ft  =  0,112  s  =  0,04 

5  5 

Subsequently,  from  Theorem  3  and  solving  the  LMI,  (40),  (41)  and  (42),  then  we  calculate  the  nonlinear 
observer  gain  matrices  L  and  ^ . 

In  this  observer,  we  apply  a  disturbance  in  the  power  supply  u  in  order  to  see  these  performances. 

H  =  10/  as  =°-°5  Q  =  5I5  Bd  =  20 

5  5  5 


"aooo5 

0 

0.0143 

0 

0 

'-  0.0574 

0 

-213.5922 

0 

0 

0.0005 

0 

0.0143 

0 

0 

-  0.0574 

0 

-213.5922 

0.0143 

0. 

0.4552 

0 

0 

G  = 

-1.7151 

0 

L  = 

2.9415 

0 

0 

0.0143 

0. 

0.4552 

0. 

0 

-1.7151 

0 

2.9415 

0 

0 

0 

0 

1.692' 

0 

0 

0 

0 

The  second  step  of  simulation  consists  of  injecting  the  obtained  numerical  values  of  the  gain  matrix 
and  the  vectors  in  an  S-function-based  Matlab  program  that  interacts  with  the  Matlab  Simulink  software  to 

simulate  the  nonlinear  system  and  the  nonlinear  00  observer. 

H 

The  simulation  results  of  the  designed  and  nonlinear  “  observer  are  presented  in  the  following. 
Figure  5  and  Figure  6  show  the  measured  and  estimated  stator  current  and  rotor  flux  components 
respectively.  The  Figure  7  and  Figure  8  show  the  measured  and  estimated  rotor  angular  velocity  and  the 
corresponding  load  torque  respectively  with  corresponding  estimation  error.  One  can  see  that  the  estimated 
state  variables  of  the  machine  follow  the  desired  trajectories. 
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Figure  5.  Measured  (red  line)  and  observed  (blue 
line)  r/-stator  current  components  and  the 
corresponding  estimation  error 


Figure  6.  Measured  (red  line)  and  observed  (blue 
line)  c/-rotor  flux  components  and  the  corresponding 
estimation  error 


The  simulation  results  of  the  designed  and  nonlinear  Hx  observer  are  presented  in  the  following. 
Figure  5  and  Figure  6  show  the  measured  and  estimated  stator  current  and  rotor  flux  components 
respectively.  The  Figure  7  and  Figure  8  show  the  measured  and  estimated  rotor  angular  velocity  and  the 
corresponding  load  torque  respectively  with  corresponding  estimation  error.  One  can  see  that  the  estimated 
state  variables  of  the  machine  follow  the  desired  trajectories. 
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Figure  7.  Measured  (blue  line)  and  observed(red 
line)  rotor  angular  velocity  of  the  (IM)  and  the 
corresponding  estimation  error 


Figure  8.  Measured  (red  line)  and  estimated  (blue 
line)  electromechanical  torque  of  the  (IM)  and  the 
corresponding  estimation  error 


According  to  the  Figure  7,  we  note  a  good  trajectory  tracking  of  the  speed  in  the  presence  of 

disturbances  and  show  the  high  performance  of  the  observer  “ .  The  simulation  results  demonstrate  that, 

despite  the  partly  unknown  transition  probabilities,  the  designed  “filters  are  feasible  and  effective, 
ensuring  the  error  systems  are  stochastically  stable. 


6.  CONCLUSION 

This  paper  gives  a  comparative  study  between  the  nonlinear  observer  based  on  circle  criterion  and 

T_T 

nonlinear  “  observer. 

The  performance  of  the  nonlinear  observers  is  evaluated  in  terms  of  their  ability  to  cope  with  model 
imperfections  and  process  uncertainties  such  as  measurement  errors  and  uncertain  initial  conditions. 

The  advantages  of  the  circle-criterion  approach  are  the  global  Lipschitz  restrictions  removing  and 
high  gain  avoiding.  However  it  introduces  linear  matrix  inequality  (LMI)  conditions. 

From  the  simulation  results,  the  proposed  “observer  has  proved  to  be  more  robust  than  the  observer  based 
circle  criterion  when  load  variations  of  the  IM  occur,  in  the  presence  of  disturbances. 

t_t 

On  the  other  hand,  the  results  obtained  using  the  nonlinear  “  observer  show  that  three 
characteristics  can  be  obtained  simultaneously.  Asymptotic  stability,  robustness  against  nonlinear  uncertainty 
and  minimized  guaranteed  H “  cost. 
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